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1 Introduction 



Many times as researchers in Combinatorics we face the task of computing 
an evaluation of the Tutte polynomial of a family of graphs or matroids. 
Sometimes this is not an easy task, or at least time consuming. Later, not 
surprisingly, we find out that a formula was known for a class of graphs or 
matroids that contains our family. Here we survey some of the best known 
formulas for some interesting families of graphs and matroids. Our hope is 
for researchers to have a place to look for a Tutte polynomial before engaging 
in the search for the Tutte polynomial formula for the considered family. 

We present along with the formulas, some explanation of the techniques 
used to compute them. This may also provide tools for computing the Tutte 
polynomials of new families of graphs or matroids. This survey can also be 
considered a companion to the book chapter Graph polynomials and their 
applications I: the Tutte polynomial by J. Ellis-Monaghan and C. Merino [25]. 
There, the authors give an introduction of the Tutte polynomial for a general 
audience of scientists, pointing out relevant relations between different areas 
of knowledge. But very few explicit calculations are made. Here we consider 
the practical side of computing the Tutte polynomial. 

However, we are not presenting evaluations, that is an immense area of 
research for the Tutte polynomial, nor analysing the complexity of computing 
the invariant for the different families. For the former we strongly recommend 
the book of D.J. A. Welsh [47], for the latter we recommend S.D. Noble's book 
chapter [36]. There are many sources for the theory behind this important 
invariant. The most useful is definitively J. Oxley and T. Brylawski book 
chapter [16]. We already mentioned [25], which also surveys a variety of 
information about the Tutte polynomial of a graph, some of it new. 

We also do not address the closely related problem of characterizing fam- 
ilies of matroids by their Tutte polynomial, a problem which is generally 
known as Tutte uniqueness. For this there are also several articles, for ex- 
ample [9, 23, 38]. 

We assume knowledge of graph theory as in [24] and matroid theory as 
in [37]. Further details of many of the concepts treated here can be found in 
Welsh [47] and Oxley and Brylawski [16]. 



2 Definitions 

Some of the richness of the Tutte polynomial is due to its numerous equivalent 
definitions; which is probably inherited from the vast number of equivalent 
definitions of the concept of matroid. In this chapter we revise three defini- 
tions and we put them to practice by computing the Tutte of some families 
of matroids, in particular uniform matroids. 

2.1 The rank- nullity generating function definition 

One of the simplest definitions, which is often the easiest way to prove prop- 
erties of the Tutte polynomial, uses the notion of rank. 

If M = {E, r) is a matroid, where r is the rank-function of M, and A O E, 
we denote r{E) — r{A) by z{A) and \A\ — r{A) by n{A), the latter is called 
the nullity of A. 

Definition 2.1. The Tutte polynomial of M, TM{x,y), is defined as follows: 
TM{x,y) = J2ix-ir^^\y-ir^^K (1) 
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Duality 

Recall that if M = {E, r) is a matroid, then M* = {E, r*) is its dual matroid, 
where r*{A) = \A\ - r{E) + r{E \ A). Because zm*{^) = nM{E \ A) and 
nM*{A) = zm{E\A), you get the following equality 

TMix,y) =TM-iy,x). (2) 

This gives the first technique to compute a Tutte polynomial. If you know 
the Tutte polynomial of M then you know the Tutte polynomial of the dual 
matroid M*. 



Uniform matroids 

Our first example is the family of uniform matroids f/r,n) where < r < n. 
Here the Tutte polynomial can be computed easily using (1) because all 
subsets of size k < r are independent, so n(-) is zero; all subsets of size k > r 



are spanning, so z{-) is zero and if a subset is independent and spanning, 
then it is a basis of the matroid. 

Tu.Jx,y) = 'f(';){x-ir-+h+± (")fe-l)" (3) 

i^Q \ / \ / i=r+l ^ ^ 

Thus, for f/2,5 we get 

TuU^,y) = {x-iy + 5{x-l) + l0 + 10{y-l) + 5{y-iy + {y-lf 
= x^ + 3x + 3y + 2y^ + y^. 

As Us^^ = (f/2,5)* 5 we get by using equation (2) 

Tc;3 ^(x, y)=x^ + 2x^ + 3x + 3y + y^. 

Matroid relaxation 

Given a matroid M = {E, r) with a subset X C E that is both a circuit 
and a hyperplane, we can define a new matroid M' as the matroid with basis 
B{M') = B{M) U {X}. That M' is indeed a matroid is easy to check. For 
example, F^" is the unique relaxation of Fj. The Tutte polynomial of M' can 
be computed easily from the Tutte polynomial of M by using equation (1). 

TM'{x,y) = TM{x,y) -xy + x + y. (4) 

Matroid relaxation will be used extensively in the last section. 

Sparse paving matroids 

Our second example can be considered a generalization of uniform matroids 
but it is a much larger and richer family. A paving matroid M = {E, r) is 
a matroid whose circuits all have size at least r. Uniform matroids are an 
example of paving matroids. Paving matroids are closed under minors and 
the set of excluded minors for the class consists of the matroid f/2,2 © f/o,i) 
see for example [33]. The interest about paving matroids goes back to 1976 
when Dominic Welsh asked if most matroids were paving, see [37]. 

Sparse paving matroids were introduced in [28, 32]. A rank-r matroid M 
is sparse paving if M is paving and for every pair of circuits Ci and C2 of 
size r we have |Ci A C2I > 2. For example, all uniform matroids are sparse 
paving matroids. 



There is a simple characterization of paving matroids which are sparse in 
terms of the sizes of their hyperplanes. For a proof, see [33]. 

Theorem 2.2. Let M be a paving matroid of rank r > 1. Then M is sparse 
paving if and only if all the hyperplanes of M have size r or r — 1. 

Note that we can say a httle more, any circuit of size r is a hyperplane. 
Conversely, any proper subset of a hyperplane of size r is independent and so 
such a hyperplane must be a circuit. Thus, the circuits of size r are precisely 
the hyperplanes of size r. 

Many invariants that are usually difficult to compute for a general matroid 
are easy for sparse paving matroids. For example, observe that if M is sparse 
paving, all subsets of size k < r are independent, and all subsets of size k > r 
are spanning. On the other hand the subsets of size r are either bases or 
circuit-hyperplanes. Thus, the Tutte polynomial of a rank-r sparse matroid 
M with n elements and A circuit-hyperplanes is given by 

TM{x,y) = j2(f)i^-^y''+Cl)+>^i^y-^-y)+i2 ifl^y-^y'^- (^) 

i^Q \ / V / i=r+l ^ ^ 

Clearly, given a r-rank paving matroid M with n elements, we can obtain 
the uniform matroid Ur^n by a sequence of relaxations from M. If M has A 
circuits-hyperplanes, equation (4) also gives (5). 

Free extension 

Another easy formula that we can obtained from the above definition involves 
the free extension, M + e, of a matroid M = {E, r) by an element e ^ E, 
which consists of adding the element e to M as independently as possible 
without increasing the rank. Equivalently, the rank function of M + e is 
given by the following equations: for X a subset of E, 

rM+e{X) =rM{X), 

and 

|r(M), otherwise. 



Again, the Tutte polynomial of M + e can be computed easily from the Tutte 
polynomial of M by using equation (1). 

TM+e{x, y) = -Tm{x, y)+ [y ) Tm(1, y). (6) 

X — 1 \ X — 1 J 

Here the trick is to notice that Tm(1,1/) = "^xiv ~ 1)\^\~^^^\ where the 
sum is over all subsets of E with r{X) = r{E), i.e., spanning subsets of M. 
The presentation given here is from [10], but the formula can also be found 
in [15]. 

For example, the graphic matroid M(i^4 \ e), has as free extention the 
matroid Qq. Because M(_ft'4 \ e) is sparse paving, by using (5), we can 
compute its The Tutte polynomial and obtain 

TK4\e{x, y) = x^ + 2x^ + X + 2xy + y + y'^. 
Now, by using (6), we get the Tutte polynomial of the free extension. 

Tgg =x^ + 3x'^ + 4x + 2xy + 4y + 3y'^ + y^. 
This result can be check by using (5), as Qq is also sparse paving. 

2.2 Deletion and contraction 

In the second (equivalent) definition of the Tutte polynomial we use a linear 
recursion relation given by deleting and contracting elements that are neither 
loops nor coloops. This is by far the most used method to compute the Tutte 
polynomial. 

Definition 2.3. If M is a matroid, and e is an element that is neither a 
loop nor a coloop, then 

Tm{x, y) = TM\e{x, y) + TM/e{x, y). (7) 

If e is a coloop, then 

Tm{x, y) = xTM\e{x, y). (8) 

And if e is a loop, then 

Tuix^y) =yTM/eix,y). (9) 
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The proof that Definition 2.1 and 2.3 are equivalent can be found in [16]. 
From this it is clear that you just need the Tutte polynomial of the matroid 
without loops and coloops. Also, if you know the Tutte polynomial of M \ e 
and M/e, then you know the Tutte polynomial of the matroid M . This way 
of computing the Tutte polynomial leads naturally to linear recursions. We 
present two examples where these linear recursions give formulas. 

The cycle C„ 

As a first example of this subsection we consider the graphic matroid M{Cn)- 
Here, by deleting an edge e from C„ we obtained a path whose Tutte poly- 
nomial is x"'~^, while if we contract e we get a smaller cycle C„_i. Thus, as 
the Tutte polynomial of the 2-cycle is x + y, we obtain 

n-l 

TcA^,y) = Y,x^ + y. (10) 

1=1 

By duality, the Tutte polynomial of the graph C* with two vertices and 
n parallel edges between them is 

n— 1 

Tc-^{x,y) = Y,y' + x- (11) 

Parallel and series classes 

As an almost trivial application of deletion-contraction we look into the com- 
mon case when you have a graph or a matroid with parallel elements. Let us 
define a parallel class in a matroid M as maximal subset X of E{M) such 
that any two distinct members of X are parallel and no member of X is a 
loop , see [37]. Then the following result is well known and has been found 
many times, see [22, 26]. 

Lemma 2.4. Let X be a parallel class of a matroid M with \X\ = p + 1. If 
X is not a cocircuit, then 

TM{x,y) = TM\x{x,y) + {y^ + y"'^ + ... + l)TM/x{x,y). (12) 

The proof is by induction on p with the case p = being (7). The 
rest of the proof follows easily from the fact that each loop introduces a 
multiplicative factor of y in the Tutte polynomial. 

A series class in M is just a parallel class in M* and by duality we get 



Figure 1: The grid graph L 



m,n- 



Lemma 2.5. Let X he an series class of a matroid M with \X\ 
X is not a circuit, then 



TM{x,y) 



(x"- 



+ xP ^ + ... + 1)Tm\x{x, y) + Tm/x{x, y) 



p + l. If 



(13) 



The rectangular lattice Lm,n 

The rectangular lattice, that we will define in a moment, is our first exam- 
ple where computing its Tutte polynomial is really far from trivial and no 
complete answer is known. The interest resides probably in the importance 
of computing the Potts partition function, which is equivalent to the Tutte 
polynomial, of the square lattice. However, even a complete resolution of 
this problem will be just a small step towards the resolution of the really 
important problem of computing the Tutte polynomial of the cubic lattice. 

Let m and n be integers, m,n > 1. The grid or rectangular lattice Lm,n 
is a connected planar graph with mn vertices, such that its faces are squares, 
except one face that it is a polygon with 2(m + n) — 4 edges. The grid graph 
Lm,n can be represented as in Figure 1 

The vertices of L^ „ are denoted as ordered pairs, such that the vertex in 
the intersection of the row i and the column j is denoted by ij. 

As example, we show a recursive formula for the Tutte polynomials of 
grid graphs -Lm.n when m = 2. For m > 3, the formulas are very complicated 
to use in practical calculations. 

The graphs ^2,15 -^2,2 and ^2,3 are shown in the Figure 2, with the labels 
that correspond to their vertices. 

For every positive integer n. The graph Q2,n is defined from the graph 



-^2,n by contraction of the edge (11, 21) of L2,„. In particular, the graph Q 
is one isolated vertex, see Figure 2 



2,1 





Figure 2: The first three terms are from the sequence {-^2,n} and the last 
three are from {Q2,n}- 

The initial conditions to construct the recurrence relations for the Tutte 
polynomial of the grid graphs L2,n are: 

i- Tq^A^.y) = 1 
ii. TL2^,{x,y) =x 
iii. Tl2_2 {x,y) = x^ + x^ + x + y 

The third condition is true because the graph L2,2 is isomorphic to the 
cycle C4, so we can use Equation (10). The formula for the recurrence relation 
to calculate the Tutte polynomial of the graphs L2,„ is: 

Tl^A^^v) = {x^ + x + l)TL2^„_,{x,y)+yTQ^^„_^{x,y) (14) 

TQ2A^^y) = ix + '^)TL2,„^i{x,y)+yTQ^^_^{x,y). (15) 

From this, we get a linear recurrence of order two. 



Tr 



{x^ + X + 1 + y)TL2,„_, + x^yTi^ 



0, 



that is easy to solve. Thus, the general formula is 

Tb^Ax^y) 



aiXi + a2 2 
"^1 



Ai — A2 
for n>2, where ai = y + x + x"^ + x^ ^ a2 



ai\2 + 02 ,„_2 



1 



A2 — Ai 
yx^ and Ai and A 2 are 

2n2\1/2' 



(16) 



{l+y + x + x'^)± ((y^ + 2y{l + x - x^) + {1 + x + x^f) 



A recursive family of graphs Gm is a sequence of graphs with the property 
that the Tutte polynomials TG^{x,y) satisfy a linear homogeneous recursion 
relation in which the coefficients are polynomials in x and y with integral 
coefficients, independent of m, see [5]. So, the sequence {-^2,n} is such a 
family and in general, {Lk^n}, for a fixed k, is a recursive family. 
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Graphic and representable matroids 

Probably the most common class of matroids, when evaluations of the Tutte 
polynomial are considered, are graphic matroids. If M is the graphic matroid 
of a graph G, then deleting an element e from M amounts to deleting the 
corresponding edge from G, and contracting e to contracting the edge. This 
is easy to do in a small graph and by using Definition 2.3 you can compute 
the Tutte polynomial. An example is shown in Figure 3 




Figure 3: An example of computing the Tutte polynomial of a graphic ma- 
troid using deletion and contraction 



As we mentioned computing the Tutte polynomial is not in general com- 
putationally tractable. However, for graphic matroids there are some re- 
sources to compute it for reasonably sized graphs of about 100 edges. These 
include Sekine, Imai, and Tani [40], which provides an algorithm to im- 
plement the recursive Definition 2.3. Common computer algebra systems 
such as Maple and Mathematica will compute the Tutte polynomial for very 
small graphs, and there are also some implementations freely available on 
the Web, such as http://ada.fciencias.unam.mx/~rconde/tulic/ by R. 
Conde or http://homepages.mcs.vuw.ac.nz/~djp/tutte/ by G. Haggard 
and D. Pearce. 

Similarly, given a matroid M = {E,r), representable over a field F, you 
can compute its Tutte polynomial using deletion and contraction. For any 
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element e & E, the matroids M \e and M/e are easily computed from the 
representation of M, see [37] . This can be automatized and there are already 
computer programs that compute the Tutte polynomial of a representable 
matroid. For our calculations we have used the one given by Michael Barany, 
for more information about this program and how to use it, see [2]. Some of 
the calculations made in the Section 4 were computed using this program. 

2.3 Internal and external activity 

Being a polynomial, it is natural to ask for the coefficients of the Tutte 
polynomial. Tutte's original definition gives its homonymous polynomial 
in terms of its coefficients by means of a combinatorial interpretation, but 
before we give the third definition of the Tutte polynomial, we introduce the 
relevant notions. 

Let us fix an ordering -< on the elements of M, say E = {ci, . . . , Cm}, 
where Cj -< Cj Hi < j. Given a fixed basis S, an element e is called internally 
active if e G S* and it is the smallest edge with respect to -< in the only 
cocircuit disjoint from S \ {e}. Dually, an element / is externally active if 
f ^ S and it is the smallest element in the only circuit contained in S'U {/}. 
We define tij to be the number of bases with i internally active elements and 
j externally active elements. In [44, 45, 46] Tutte defined Tm using these 
concepts. A proof of the equivalence with Definition 2.1 can be found in [6]. 

Definition 2.6. If M = {E,r) is a matroid with a total order on its ground 
set, then 

TMix,y) = ^tijx'y^ . (17) 

In particular, the coefficients tij are independent of the total order used on 
the ground set. 

Uniform matroids again 

Our first use of (17) is to get a slightly different expression for the Tutte 
polynomial of uniform matroids. This time we get 



-..(-)^|:C;!;>'-i:(:;:::;)^ 
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when < r < n, while Tu^^^{x,y) = x" and Tug^{x,y) = y". This can also 
be established by expanding (3). 



Paving matroids 

Paving matroid were defined above. The importance of paving matroids is 
its abundance, meaning, most matroids of up to 9 elements are paving. This 
was checked in [32] and it has been conjecture in [31] that this is indeed true 
for all matroids, that is, when n is large, the probability that you choose 
a paving matroid uniformly at random among all matroids with up to n 
elements is approaching 1. 

Now, in order to get a formula for the Tutte polynomial of paving ma- 
troids, we need the following definition from [37]. Given integers k > 1 and 
m > 0, a collection T = {Ti, . . . ,Tk} of subsets of a set E, such that each 
member of T has at least m elements and each m-element subset of E is con- 
tained in a unique member of T, is called an m-partition of E. The elements 
of the partition are called blocks. The following proposition is also from [37]. 

Proposition 2.7. // T is an m-partition of E, then T is the set of hyper- 
planes of a paving matroid of rank m + l on E. Moreover, for r > 2, the set 
of hyperplanes of every rank-r paving matroid on E is an (r — 1) -partition of 
E. 

Brylawski gives the following proposition in [13], but we mentioned that 
his proof does not use activities. 

Proposition 2.8. Let M be a rank-r matroid with n elements and Tutte 
polynomial '^ijtijX^y^. Then, M is paving iff tij = for all {i,j) > (2,1). 
In addition if the (r — l)-partition of E has h^ blocks of cardinality k, for 
k = r — 1, . . .n, then 

I Tl — t — 1 \ 

tio = { . ] for all i > 2; 

y r — t J 
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and for all j > 0, 



'li 



toj 



X] ( ^ o )^fc+i+r-i, and 



fc=0 



r-2 



1 )2^[ r-l l^'^+J+'-i- 



r — 



Note that when M is the uniform matroid f/r,n, M is a paving matroid 
with (r — l)-partition all subsets of size r — 1. Thus the above formula give 
us equation (18). 




2 4 5 

Figure 4: The matroid Rg 



As a example consider the matroid -Rg with geometric representation given 
in Figure 4. This matroid is paving but not sparse paving. The blocks 
correspond to lines in the representation. There are 7 blocks of size 3, 2 
blocks of size 4 and 3 blocks of size 2, corresponding to trivial lines that 
do nor appear in Figure 4. Using the above formulas we get that the Tutte 
polynomial of Rg is 



TR,ix,y) 



X 



6x^ + 8x + llxy + 2xy^ + Sy + 13y^ + lOy^ + 6y^ + 3y^ + y^ 



Catalan matroids 

Our final example in this subsection is a family of matroids that although 
simple, it has a surprisingly natural combinatorial definition together with a 
nice interpretation of the internal and external activity. 
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Consider an alphabet constituted by the letters {A^, E}. The length of a 
word w = WiW2---Wn is n, the number of letters in w. Every word can be 
associated with a path on the plane Z^, with a initial point A = (mi, ni) and 
a final point B = (m2,n2), mi < 1712 and rii < n2. The letter N is identified 
with a north step and the letter E with a east step, such that the first step 
oi A = (0, 0) to B will be A^ = (0, 1) or E = (1, 0). 

Let A and 5 be a couple of fixed points oil?. A lattice path is a path in 
1? from A to -B using only steps N or E. The lattice paths from Ato B have 
m steps ii^ and n steps A^. Let P and Q be two lattice paths, with initial 
points (xo. Up) and (xq, t/g) of P and Q respectively. If yp < yq for every xq 
in [mi, 7722]; then the set of lattice paths in this work are lattice paths from A 
to B bounded between P and Q, that we called PQ-bounded lattice paths. 

The lattice paths P and Q can be substituted by lines. If P is the line 
y = 0, Q the line y = x and m = n, then the number of PQ-bounded lattice 
paths from (0, 0) to (n, n) is the n-th Catalan number: 

1 / 2n 



n+l 



n 



Denote by [n] the set {1, 2, ..., n}. Consider the word w=Wi, W2, ..., Wm+n 
that corresponds to a PQ-bounded lattice path. We can associate to w a 
subset Xyu of [m+n] by defining that i G X^ iiwi = N. X^ is the support set 
of the steps N in w. The family of support sets X^ of PQ-bounded lattice 
paths are the set of bases of a transversal matroid, denoted M[P, Q] with 
ground set [m + n], see [10]. 

If the bounds of the lattice paths are the lines y = and y = x, the 
matroid is denoted by M„ and is named Catalan matroid. Note that M„ has 
a loop that corresponds to the label 1 and a coloop that corresponds to the 
label 2n, and that M„ is autodual. 

The fundamental result to compute the Tutte polynomial of Catalan ma- 
troids is the following, see [10], where here, for a basis B, we denote by i{B) 
the number of internally active elements and by e{B) the externally active 
elements. 

Proposition 2.9. Let B E B be a basis of M[P,Q] and let wb be the PQ- 
bounded lattice path associated with B. Then i{B) is the number of times wp 
meets the upper path Q in a north step and e{B) is the number of times wp 
meets the lower path P in an east step. 
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Figure 5: Bases of M3 with their internal and external activity. 



Then the Tutte polynomial of the Catalan matroids M„ for n > 1, is: 



Tm„{x 



,y)=j: 



i + j 



i,j>0 



n — 1 



2n — t — ] — I ■ ,1 

,1 x^y ■ 



(19) 



Note that the coefficient of x'^y^ in the Tutte polynomial of the matroid 
M„ depends only on n and the sum i+j. 

In the Figure 5 are shown the bases of M3 with its internally and external 
activities. Thus the Tutte polynomial of M3 is 

Tms {x, y) = x^y + x^y + x^y^ + xy^ + xy^. 



3 Techniques 



We move now to present three techniques that are sometimes useful to com- 
pute the Tutte polynomial, when the initial trial with the above definitions 
are not successful. 
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3.1 Equivalent polynomials: The coboundary polyno- 
mial 

For some families of matroids and graphs it is easier to compute polynomials 
that are equivalent to the Tutte polynomial. One such polynomial is Crapo's 
coboundary polynomial, see also [14]. 

XM{X,t)= J2 ^""xm/x(A) (20) 

where J-'m is the set of flats of M and Xn{^) is the characteristic polynomial 
of the matroid A^. The characteristic polynomial of a matroid M is defined 
by 

X^(A) = E(-l)""^^^^^- (21) 

ACE 

The Tutte and the coboundary polynomial are related as follows 

TM{x,y) = ———^ XM{{x-l){y-l),y). (22) 

XM(A,t) = (t - l)'-(^) TMi ^^^_~\ t). (23) 

When M is a graphic matroid M{G) of a connected graph G, the charac- 
teristic polynomial is equivalent to the chromatic polynomial, that we denote 
Xg- 

XgW = Axm(g)(A). 

And also, in this case, the coboundary is the bad-colouring polynomial. 
The bad colouring polynomial is the generating function 

5G(A,t) = 5^6,(G;A)t^ 

3 

where bj{G; A) is the number of A-colourings of G with exactly j bad edges. 
Note that when we set t = we obtain the chromatic polynomial. Then, the 
polynomials are related as follows: 

B{G;X,t) = XxM(G){\t). (24) 
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The g-cone 

We base this section on the work of J. E. Bonin and H. Qin, see [11]. 

Definition 3.1. Let M be a rank-r simple matroid representable over GF(q). 
A matroid N is a q-cone of M with base S and apex a if 

1. the restriction PG{r-^q)\S of PG{r,q) to the subset S is isomorphic to 
M, 

2. The point a G PG{r, q) is not contained in the closure of S, GLpG(r,q){S), 
in PG{r, q), and 

3. The matroid N is the restriction of PG{r,q) to the set[J^^gGlpG(r,q)i{0',s}). 

That is, one represents M as a set S of points in PG{r, q) and constructs 
A^ by restricting PG{r, q) to the set of points on the hnes joining the points 
of S to the fixed point a outside the hyperplane of PG{r, q) spanned by S. 
The basic result is by Kung who proved in [29]. 

Theorem 3.2. For every q-cone N of a rank-r simple matroid M, we have 

X;v(A) = (A-l)g^XM(A/g). 

To extend this result to x the authors in [11] identify the flats F of N and 
the contractions N/F of A^ by these fiats. They prove the following formula 

Xn{\ t) = txMiK f^) + q''{X - l)XM(A/g, t). (25) 

By using equation (22) we get a formula for the Tutte polynomial of A^ in 
terms of the Tutte polynomial of M. 

Theorem 3.3. If M is a rank-r matroid representable over GF{q) and N is 
a q-cone of M, then 

rp f X yJy'^-'^Y rp f ix-l){y-l) q''{xy-x-y) x + 1 
TN{x,y) = 7 TvTI^^ ~a — ^ +l>r + -. Tm{ +1,1/ • 

Thus, for example, the 2-cone of the 3-point line is the matroid PG{2, 2). 
The Tutte polynomial of the former matroid is just x'^ -\-x -\-y, then by using 
the above formula we get the Tutte polynomial of PG{2, 2) to be 

x^ + 4x^ + 3a; + 7xy + By + Qy^ + 3y^ + y*. 

As the g-cone of PG{r — 1, q) is PG{r, q), this method can be used to com- 
pute the Tutte polynomial of any PG{r,q), but we do this later using the 
coboundary polynomial directly. 
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Complete graphs 

Given the apparent simplicity of many of the formulas for invariants in com- 
plete graphs, like the number of spanning trees or acyclic orientations, it 
is not so straightforward to compute the whole Tutte polynomial of com- 
plete graphs; many researches, however, have tried with different amounts of 
success. The amount of frustration after failing to compute this seemingly 
innocuous invariant of many of our colleagues was our original motivation 
for writing this survey paper. 

Using the exponential formula, see Stanley [42], we can give an exponential 
generating function for the Tutte polynomial of the complete graphs. Let us 
denote the vertices of Kn by V and by -B„ its bad-colouring. To compute Bn, 
observe that any A-colouring partitions the vertices V into A color classes of 
subsets Vi of vertices each of cardinality nj, for i = 1 . . .\. So, we have that 
ni + . . . + n\ = n. The number of bad edges with both ends in the set Vi is 
t^ 2V. Thus, by the exponential formula we get the following formula 




= l + ^i?„(A,t)-. (26) 

n>l 

Now, we can use equation (24) to get 

E'<">;^)' = l + ^EX»(A.05^ (27) 

K.n&i J n>l 

Let Tn{x,y) be the Tutte polynomial of Kn- Tutte in [46] and Welsh 
in [48] give the following exponential generating function for Tn{x,y) that 
follows from the previous equation and equation (22). 

„\ (-i)(J/-i) 

E^^^^^ =l + (x-l)5^(y-l)'^T„(a:,i/)^. (28) 

■^ — ^ n! / ^ — ^ n! 

\n>0 / n>l 

Even that the previous formulas seem difficult to handle, equation (27) 
it is quite easy to use in Maple or Mathematica. For example, the following 
program in Maple can compute T3o(x, y) in no time. 

Coboundary : = proc(n,q,v) local i,x,g; 

g:=x->(add(v"(i*(i-l)/2)*x"i/i! ,i=0. .n))"q; 



simplify (eval(diff(g(x) ,x\$n) ,x=0)/q) ;end proc; 
T:= proc(n,x,y) 

simplifyC (l/((y-l)-(n-l)))*subs({q=(x-l)*(y-l) ,v=y}, 
Coboundary(n,q,v)) ) end proc; 

By computing T^{x,y) we obtain 

T5(x, y) = y^ + Ay^ + X* + 5 xy^ + 10 y'^ + 6 x^ + 10 x'^y + 15 xy'^ + 15 y^ 
+ llx'^ + 20xy+15y^ + 6x + 6y. 

Complete bipartite graphs 

It is just natural to use this technique for complete bipartite graphs Kn,m 
with similar results as above. This time let the vertex set be Vi U V2 and 
denote by Bn^m hs bad-colouring polynomial. To compute -B„,m, observe 
that any A-colouring partitioned the vertices in subsets V^ U V^ each of 
cardinality niirii. The number of the bad edges with both ends colour i is 
nifrii. Thus, by the exponential formula we get the following. 



n! ml I ^-^ n\ ml 

(n,m)7^(0,0) (29) 

= 1 + A > Xn,m{\t) — —r. 

^-^ n! ml 

{n,m)^(0,0) 

Thus, a formula for the Tutte polynomial of the bipartite complete graph 
can be found in a similar way as before. Let Tn,m{x, y) be the Tutte polyno- 
mial of Kn,m- The following formula is from Stanley's book [42], see also [30] 




•^— ^ n I ml / 



yji ^m 



nl ml 

^(n,m)6N2 , 

^-^ nl ml 

(n,m)eN2 
(n,m)^(0,0) 
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Figure 6: The proyective plane PG{2,3) with a marked hyperplane 



As before, equation (29) is quite easy to use in Maple to compute Tn^mix, y) 
for small values of n and m. In this way we get the following. 



T3_3(x,y) = x^ + 4 x^ + 10x^ + 9x^1/+ llx^ + 6x|/^ 
+ 15 xy + 5 X + y^ + 5 y^ + 9 y^ + 5 y. 



(31) 



Projective geometries and afRne geometries 

The role played by complete graphs in graphic matroids is taken by projective 
geometries in representable matroids. Herein lies the importance of projective 
geometries. Even though a formula for their Tutte polynomial has been 
known for a while, and has been discovered at least twice, not much work 
has been done in the actual combinatorial interpretations for the value of the 
different evaluations of the Tutte polynomial. 

For this part we follow Mphako [34], see also [3]. For all non-negative 
integers m and k we define the Gaussian coefficients as 



m 
k 



k-U 



(g'"-l)(g™-g)---(g"^-g 
{q'' — l){q'' — q) ■ ■ ■ {q'^ — q^^^) 



(32) 



Note that [™] =1 since by convention an empty product is 1 and [°] =0. 

Let PG{r — l,q) be the (r — 1) dimensional projective geometry over 

GF{q). Then, as a matroid, it has rank r and [^] elements. Also, every 

rank-A; flat X is isomorphic to PG{k — 1, q) and the simplification of M/X 
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is isomorphic to PG{r — k — l,q). The number of rank-A; flats is [^] . The 
characteristic polynomial of PG{r — l,g) is known to be, see [16], 



r-l 



XPG(r-i,g)(A) = JJ(A-g* 
Thus, using equation (20) we get 



i=0 



XPG(r-l,,)(A,t) = 5^t['] 



r— fc— 1 



A;=0 



k 



n (^-^^)- 



(33) 



j=0 



The matroid PG{2, 3), has a geometric representation given in Figure 6, 
and it is isomorphic to the unique Steiner system S'(2,4, 13). The matroid 
is paving but not sparse paving and has the following representation over 
GF{3). 

1 2 4| 3 5 6 7 8 9 10 11 12 13 
1 0| 1 2 2 1 1 1 1 2 
01 01 1110111 1 1 



11 1 2 1 2 1 







To compute its Tutte polynomial we could use the program in [2] or use the 
above formula to get 

Xpg{2,3)(A, t) = (A - 1)(A - 3)(A - 9) + 13t(A - 1)(A - 3) + 13t^(A - 1) + 1'^ 

Finally the Tutte polynomial of PG{2, 3) is obtained from the coboundary 
by a substitution X = {x — l){y — 1) and t = y and by multiplying by the 
factor l/{y — 1)^ 



T, 



PG{2,3) 



'x,y) 



x^ + lOa;^ + 13xy^ + 2Qxy + IQx + IQy + 32y^ + 36y^ 
+28y^ + 21y^ + I5y^ + lOy'^ + 6y^ + 3y^ + y^°. (34) 



Every time we have a projective geometry PG{r—l, q) we can get an affine 
geometry^ AG{r — 1, g), simply by deleting all the points in a hyperplane of 
PG{r — l,g). For example, by deleting all the red points from PG'(2,3) in 
Figure 6 we get AG{2, 3) with geometric representation in Figure 7. 

The situation to compute the Tutte polynomial is very similar for the 
r-dimensional affine geometry over GF{q), AG{r,q). It has rank r + 1 and 
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Figure 7: Afine plane AG{2, 3) 



q"^ elements. Any flat of X rank k + 1 is isomorphic to AG{k,q) and the 
simplification of AG{r, q)/X is isomorphic to PG{r — k — l,q). The number 
of rank-A; flats is g*""^ [^] . The characteristic polynomial of AG{r,q) is 
known to be, see [16], 



fc-i 



XAG(.,.)(A) = (A-l)J](-l)^y-^n(^^ 



fc=0 



1=0 



Now, using the result of Mphako we can compute the coboundary polynomial 
of ^^(r,^). 



XAG(r,g)(A,t) = XAG{r,q)W + J^^^'? 



r—k 



k=0 



r—k—1 



n (^-^^)- 



(35) 



j=0 



The afiine plane AG{2, 3) is isomorphic to the unique Steiner triple system 
5'(2,3, 9), as every line contains exactly 3 points and every pair of points 
is in exactly one line. S{2, 3, 9) is a sparse paving matroid with 12 circuit- 
hyperplanes so by using equation (5) we could compute its Tutte polynomial. 
However, we use the above formula to obtain the same result. 

XAG{2,3)(A, t) = (A - 1)(A2 - 8A + 16) + 9t(A - 1)(A - 3) + I2t\\ - 1) + f. 

Again, the Tutte polynomial of AG{2, 3) is obtained from the coboundary by 
a substitution A = {x — l){y — 1) and t = y and by multiplying by the factor 
1/iy-ir 



T, 



AG(2,3) 



[x, y) =x^ + 6x^ + 12xy + 9x + 9y + ISy^ + lOy^ + 6y^ + 3y^ + y^. 
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3.2 Transfer-matrix method 

Using formula (1) quickly becomes prohibitive as the number of states grow 
exponentially with the size of the matroid. You can get around this problem 
when you have a family of graphs that are constructed using a simple graph 
that you repeat in a path-like fashion; the bookkeeping of the contribution 
of each state can be done with a matrix, the update can be done by matrix 
multiplication after the graph grows a little more. This is the essence of our 
second method. 

The theoretical background of the transfer-matrix method, taken from 
[41], is described below. 

A directed graph or digraph G is a triple (V, E, 0), where V = {vi, . . . , Vp} 
is a set of vertices, E is a. finite set of directed edges or arcs, and is a map 
from E to V X V. If 0(e) = {u, v), then e is called an edge from u to v, with 
initial vertex u and final vertex v. A directed walk F in G of length n from u 
to f is a sequence ei, . . . , e„ of n edges such that the final vertex of Cj is the 
initial vertex of Cj+i, for 1 <i <n — 1. 

Now let w : -E — > -R be a weight function on E with values in some 
commutative ring R. If F = ei, . . . , e„ is a walk, then the weight of F is 
defined by w{T) = w{ei) ■ ■ ■ w(e„). For 1 < i,j < p and n G N, we define 

r 

where the sum is over all walks F in G of length n from Vi to Vj. In partic- 
ular, Aij{0) = 6ij. The fundamental problem treated by the transfer- matrix 
method is the evaluation of Aij{n). The idea is to interpret Aij{n) as an 
entry in a certain matrix. Define a p x p matrix D = (Dij) by 



D^. = E 



wie] 



where the sum is over all edges e satisfying that its initial vertex is Vi and 
its final vertex is Vj. In other words, Dij = Aij{l). The matrix D is called 
the adjacency matrix of G, with respect to the weight function w. 

Theorem 3.4. Let n G N. Then the {i,j)-entry of D^ is equal to Aij{n). 
Here we define D^ = Ip even if D is not invertible, where Ip is the identity 
matrix. 

Proof See [41]. D 
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Rectangular lattice again 

The transfer-matrix method gives us another way to compute Tl^ ^ {x, y) for 
a fixed width m at point {x,y) which is described in Calkin et. al. [17]. In 
this case we have the same restriction as before, a fixed width m for small 
values of m, but it has the advantage of being easily automatized. 

Theorem 3.5 (Calkin et. al. 2003). For indeterminates x andy and integers 
n,m > 2, m fixed, we have 

Tl^Jx + l,y + l)= x-^-'Xl . (A^)"-i ■ 1, 

where X.^, a vector of length Cm, and Am, a Cm x Cm matrix, depend on x,y 
and m but not n. And 1 is the vector of length Cm with all entries equal to 1. 

The quantity Cm is the m-th Catalan number, so the method is just 
practical for small values of m. Computing the vectors X^ and the matrix 
Km can be easily done in a computer. 

For example, for L2,„ we get that T^^^^x + 1, y + 1) equals 



x'--\x-\l)[ U 



x"^ + 3x'^ + yx~'^ 1 + 2x'^ \ (I 



X 



Wheels and whirls 



2x~'^ + a;-^ 1 + 2x-^ + x" 



■2 



The transfer matrix method takes a nice turn when combined with the physics 
idea of boundary conditions. In this case more lineal algebra is required but 
the method is still suitable to use in a computer algebra package. 

A well-know family of self-dual graphs are wheel graphs, Wn- The graph 
Wn has n + 1 vertices and 2n edges, see Figure 8. The rim of the wheel 
graph Wn is a circuit-hyperplane of the corresponding graphic matroid, the 
relaxation of this circuit-hyperplane gives the matroid W^", the whirl matroid 
on n elements. In [19], S.-C. Chang and R. Shrock using results from [39] 
compute the Tutte polynomial of Wn- 

Tw„{x,y) = ^[{l+x + y) + i{l + x + yf-AxyY^'r (36) 

+ — [(1 +x + y)-{{l + x + yy- 4xyy/^r + xy - x - y - 1. 
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Figure 8: The n-wheel 



From this expression is easy to compute an expression for the Tutte poly- 
nomial of whirls using equation (4). 

Tw<x,y) = l-[(l + x + y) + {{l+x + yf-AxyY/'r (37) 

+ ^[{l + x + y)-{{l + x + y)'-4xyy/'r-l. 

Now, the way S.-C. Chang and R. Shrock compute T]v^{x,y) is by using 
the Potts model partition function together with the transfer matrix method. 
Remarkably, the Tutte polynomial and the Potts model partition function 
are equivalent. But rather than defining the Potts model we invite the reader 
to check the surveys in [49, 4]. Here we show how to do the computation 
using the coboundary polynomial and the transfer method. 

Let us compute the bad colouring polynomial for Wn when we have 3 
colours. For this we define the 3x3 matrix D^. 



e 


1 


1 


t 


t 


1 


t 


1 


t 



The idea is that the entry ij of the matrix (-Da)"^ = -D3 will contain all 
the contributions of bad edges to the bad colouring polynomial for all the 
colourings a, with a{h) = 1, o"(l) = i and a{n) = j for the fan graph F„, 
that is obtained from Wn by deleting one edge from the rim. Then, to get the 
bad-colouring polynomial of F„ we just add all the entries in D3 and multiply 
by A = 3. To get the bad-colouring polynomial of Wn we take the trace of 
D^ and multiply by A = 3. This works because, by taking the trace you are 
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considering just colourings with the same initial and final configuration, that 
is you are placing periodic boundary conditions. 

Thus, for example, the trace in D| is t^ + 8t^ + 6t^ + 12t, so Bw^i^yt) = 
3t^ + 24t^ + 18t^ + 36t. By computing the eigenvalues of D^ we obtain the 
bad colouring polynomial of Wn with A = 3. 



-Sw'„(3,t) 



-[(t^ + t + 1) + (t^ - 2t^ - t^ + lot + 1 



.l/2in 



-[{t^ + t + 1) + (t^ - 2t^ -e + lot + If'^f + 3(t - 1)". 



For arbitrary A we need to compute the eigenvalues of the matrix of A x A 
given below 



(e 1 1 
t t 1 
t 1 t 



i\ 

1 

1 



yt 1 1 ■■■ tj 
The eigenvalues of the matrix are 



ei,2 = -[(t' + t + A-2)±(t' 



2r 



(2A - 5)t2 + (6A - 8)t + (A - 2fY/\ 



each with multiplicity 1 and 63 = (t — 1) with multiplicity A — 2. Thus the 
bad colouring polynomial of Wn is A(e" + e2 + (A — 2)63). To get formula (36) 
we just need to make the corresponding change of variables as seen at the 
beginning of the previous section. 

Other examples that were computed by using this technique include 
Mobius strips, cycle strips and homogeneous clan graphs. For the corre- 
sponding formulas you can check [39, 18, 19, 20, 21]. 

3.3 Splitting the problem: 1, 2 and 3-sunis 

Our last technique is based on the recurrent idea of splitting a problem that 
otherwise may be too big. The notion of connectedness in matroid theory 
offers a natural setting for our purposes, but we do not explain this theory 
here and we direct the reader to the book of Oxley [37] for such a subject. 
Also, for 1, 2 and 3-sums of binary matroids the book of Truemper [43] is 
the best reference. 
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However, we would like to explain a little of the relation of 1, 2 and 3- 
sums and the splitting of a matroid. For a matroid M = {E,r), a partition 
{X, Y) of E is an exact k-separation, for k a positive integer, if 

min{|X|,|r|} >k 

and 

r{X)+r{Y)-r{M) = k - I. 

Now, a matroid M can be written as a 1-sum of two of its proper minors if 
and only if M has an exact 1-separation; and M can be written as a 2-sum 
of two of its proper minors if and only if M has an exact 2-separation. The 
situation is more complicated in the case of 3-sums and we just want to point 
out that if a binary matroid M has an exact 3-separation {X,Y), with |X|, 
|y| > 4, then there are binary matroids Mi and M2 such that M = Mi©3M2. 

1-sum 

For n matroids Mi, M2, . . ., M„, on disjoint sets Ei, E2,. . ., En the direct 
sum Ml © M2 © • • • © Mn is the matroid {E,X) where E is the union of the 
ground sets and X = {/i U /2 U . . . U /„ : Jj G T{Mi) for all i in {1, 2, ... , n}}. 
Directly from (1) it follows that the Tutte polynomial of the 1-sum is 
given by 

n 

TMi(BM2®-(BM„{x,y) = YlTM,{x,y). 

i=l 

2-suin 

Let Ml = (^i,Xi) and M2 = (£^2,2^2) be matroids with £^1 n £^2 = {p}- If P 
is not a loop or an isthmus in Mi or M2, then the 2-sum Mi ©2 M2 of Mi 
and M2 is the matroid on £^1 U -£2 \ {p} whose collection of independent sets 
is {h Uh: he Xi, I2 e X2, and either h U {p} G Xi or h U {p} G X2} 

In [1], [12] and [37], we find recursive formulas for computing the Tutte 
polynomial of the matroid Mi ©2 M2, in term of the matroids Mi and M2. 
Here we present the formula given in [1] for TMi®2M2{.x,y). 



Tmi®2M2 — X/vf /r TMi\r\ 

xy — X — y ^ -' 
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X- 1 - 1 

-1 y-1 



TM2/r 
TM2\r 
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where here we omit the variables x, y of each Tutte polynomial to avoid a 
cumbersome notation. 

As an example, let us take Rq that is the 2-sum of f/2,4 with itself. Observe 
that the election of the base point is irrelevant as its automorphism group is 
the symmetric group. The matroid f/2,4/''" — f/1,3 and f/2,4 \ r = f/2,3. Then, 
by using (38) we get that Tu2 402^^2 4 equals 



[y'^ + y + X x'^ + X + y~\ 



X- 1 - 1 
-1 y-1 



y"^ + y + X 

x^ + X + y 



xy — X — y 
Thus, we obtain that T^g {x, y) = x^ + 3x^ + Ax + 2xy + % + Sy"^ + y^. 

3-sum 

The best known variant of a 3-sum of two matroids is called A-sum, see [43]. 
For 3-conected matroids Mi and M2, the A-sum is usually denoted M = 
Ml ©1^ M2. When Mi and M2 are graphic matroids with corresponding 
graphs being Gi and G2, the graph G = Gi ©^ G2 is obtained by identifying 
a triangle of Ti of Gi with a triangle T2 of G2 into a triangle T, called 
the connector triangle. Finally, G is obtained by deleting the edges in the 
connector triangle. 

Here, we present a formula to compute the Tutte polynomial of the A- 
sum of Ml and M2, in terms of the Tutte polynomial of certain minors of the 
original matroids. The expression for Tj^^^i\j^^{x,y) was taken from [1] and 
its proof is based on the concept of bipointed matroids that is an extension 
of the pointed matroid introduced by Brylawski in [12]. Also, the work in [1] 
is more general as the author gives an expresion for the Tutte polynomial of 
a certain type of general parallel connection of two matroids. 

Let Ml and M2 be two matroids defined on Ei and -^2, respectively. Let 
T be equal to Ei D E2 = {p, s, q}, a 3-circuit and A^ = Mi\T = M2IT. We 
requiere that in Mi there exist circuits f/iU{s} with f/i C Ei\T and f/2U{p} 
with U2 '^ Ei\ T; similarly, we need that in M2 there exist circuits Vi U {s} 
with ViC E2\T and V2 U {p} with \4 C E2 \ T. 

For i = 1, 2, ..., 5, we consider the following 5 minors Qi of Mi, on Ei \ T, 
and 5 minors Pi of M2, on E2 \ T. 

Qi = Mi\p\s\q,Q2 = Mi\p/s\q,Q3 = Mjp \s\q,Q4 = Mi/p/s/q 
and Q5 = Mi\p\ s/q. 
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Pi = M2\p\s\q, P2 = M2\p/s\q, P3 = M2/p\s\q, P4 = M^/p/s/q and 
P5 = M2\p\s/q. 

We take the vectors over Z[x,y]: 

1^ = [Tqi,Tq2,Tq,^,Tq^,Tq^,], 

again here we omit the variables x, y of each Tutte polynomial to avoid a 
cumbersome notation. 

Finally, the formula for the 3-sum (A-sum) of Mi and M2, it is as follows. 



T, 



Mie^Mz 



(x,|/) = (g"")* 



— 1 — X — y + xy 



C 



P 



(39) 



where the matrix , C is given by: 



— z— J/+XJ/ — x— j/+a;j/ —x~y+xy —x—y+xy 

l—y -| 1 1—x 

—x—y+xy -~x~y+xy —x—y+xy 

l—y 1 1 1— a: 

— X— jz+xj/ — z— j/+a;j/ —x—y+xy 

2 1-x 1-x (1-a:)^ 

—x—y+xy —x—y+xy —x—y+xy —x—y+xy 

l-y 1 1 1-x 

—x—y+xy —x—y+xy —x—y+xy —x—y+xy 



l-y 



-x—y+xy 

1 
-x—y+xy 

1 
-x—y+xy 

1-x 
-x—y+xy 

1 



As an example let us take F^ that is the 3-sum of Fj and Pf along a 
3-circuit. In this case we have for P7 the following table, where in the first 
column we present the 5 minors we need, then the second column has the 
corresponding matroid, and the third column has the corresponding Tutte 
polynomial of that matroid. 
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Minor 


Matroid 


polynomial 


F7\p\s\q 


f/3,4 


x^ + x^ + X + y 


Fj \p/s\q 


f/1,2 © f/1,2 


{x + yf 


Fr/p\s\q 


f/1,2 © f/1,2 


{x + yf 


Fi/p/s/q 


f/1,4 


y^ + y'^ + y + X 


Ft\p\ s/q 


f/1,2 © f/1,2 


(x + y)^ 



Similarly for F7 we have the corresponding table: 



Minor 


Matroid 


polynomial 


F,-\p\s\q 


f/3,4 


x^ + x'^ + X + y 


F,~\p/s\q 


C3 plus a paralell edge 


x'^ + X + xy + y + y'^ 


F^/p\s\q 


f/1,2 © f/1,2 


{x + yf 


F^/p/s/q 


f/1,4 


y^ + y^ + y + x 


F^\p\ s/q 


f/1,2 © f/1,2 


{x + yf 



By using (39) with the values of the above tables we get Tp^ {x, y) = 
x'^ + 4x^ + lOx^ + 8a; + Vlxy + Sy + lOy^ + Ay^ + y^. 

There is a general concept oik — sum for matroids and a formula exists 
for this general notion, however the formula is quit intricate, so we refer the 
reader to the original paper of Bonin and de Mier in [8]. 

Thickening, stretch and tensor product 

Given a matroid M and a positive integer k, the matroid M^''^ is the ma- 
troid obtained from M by replacing each non-loop element by k parallel 
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elements and replacing each loop by k loops. The matroid M^'^^ is called 
the k-thickening of M. In [16] the following formula is given for the Tutte 
polynomial of M^''^ in terms of that of M. 



Tmw {x, y) 



,fc-l I „,k-2 



(/-i + /-2 + ... + ^ + i)KM)y. ^y +y 



y + x fc, (40) 



M\ t._i h—9 I I I 1 ' ^ 

yk I J^yH ^_|___|_y_|_X 

A proof by using the recipe theorem can be read in the aforementioned 
reference. Here we hint a simple proof by noticing that any flat of M^^'> is 
the k-thickening of a flat of M . Thus, by equation (20) 

xM(^)(g,t)=xM(g,t'=). 

And thus, by using equations (23) and (22) we obtained the formula. 

The dual operation to A;-thickening is that of k- stretch that is defined 
similarly. The matroid M(fc) is the matroid obtained by replacing each non- 
isthmus of M by A; elements in series and replacing each isthmus by k isth- 
muses. The matroid M(fc) is called the k-stretch. It is not difficult to prove 
that M(fc) = ((M*)*^*'^)* and so, we obtained the corresponding formula for 



Tm, 



(fe)- 
^M(fe) {x, y) = 



x''-'^ +x^-^ + ... + x + y. (41) 



^k l_|_a;fe 2_|__^_|_2;_)_1 

More generally we have the following operation called the tensor product. 
A pointed matroid A'"^ is a matroid on a ground set which includes a distin- 
guished element, the point d, which will be assumed to be neither a loop or 
coloop. For an arbitrary matroid M and a pointed matroid A^, the tensor 
product M Ad is the matroid obtained by taking a 2-sum of M with A^^ at 
each point e of M and the distinguished point d of A^^. The Tutte polynomial 
of M ® Na, where M = {E, r) is then given by 

T (^ n,\ - f\E\-r(E) r{E)rj. A^ - 1) f + Q f + {y - l)g . , . 

iM(^NAx,y> - J 9 J-mK , 7 )■, l42j 

where / = f{x,y) and g = g{x,y) are polynomials which are determined by 
the equations 

{x-l)f{x,y)+ g{x,y) = TN^\dix,y) 
f{x,y) + {y-l)g{x,y) = TN,/d{x,y). 
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The proof of the formula uses a generahzation of the recipe theorem to 
pointed matroids and can be found in [15], here we follow the exposition 
in [16, 27]. Observe that when a matroid A^ has a transitive automorphism 
group the choice of the distinguished point d is immaterial. Thus, if N is 
Uk^k+i, k > 1, we get the /c-stretch and if N is Ui^k+i, k > 1, we get the 
/c-thickening. 

4 Aplication: Small matroids 

Let us put these techniques in practice and compute some Tutte polynomials 
for matroids with a small number of elements, these are matroids from the 
appendix in Oxley's book [37]. We will try, whenever possible, to check 
the result by using two techniques. Soon, the reader will realize that a fair 
amount of the matroids considered are sparse paving and that computing the 
Tutte polynomial for them is quite easy. 

Matroid f/2 4 

The Tutte polynomial of the uniform matroid f/2,4 can be computed using 
equation (3). 

Tu,,,{x,y)=x^ + 2x + 2y + y\ (43) 

Also this matroid is the 2-whirl W"^ so you can check this result using 
equation (37) 



Figure 9: The uniform matroid f/2,4 



Matroids W3, W^, W^ and W^ 

The 3-wheel M(W^3), which is isomorphic to the graphic matroid M(i^4), is 
a sparse paving matroid with A = 4 circuit-hyperplanes, so by using (5) we 
get 

Tw3 (x, y) =x^ + 3x^ + 2x + Axy + 2y + 3y^ + y^. (44) 
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Of course, the above polynomial can be checked using (36). The only 
relaxation of M{W3) is the 3-whirl W^, so by using (4) we obtain its Tutte 
polynomial. This can be checked by using (37). 



T\Y3{x, y) = x^ + 3x^ + 3x + 3xy + 3y + 3y^ + y^. 



(45) 



The 4-wheel, MiWi), and 4- whirl, W^, are matroids whose geometric 
representation are shown in Figures 10 and 11. Their Tutte polynomial can 
be computed using (36) and (37). 



TwA^^y) 



{l + x + y+^{l + x + yf- 


- AxyY 


16 


{l + x + y-^{l + x + yY- 


- AxyY 



16 



+ 



+ xy — X — y — 1 



x^ + Ax^ + Qx^ + 3x + Ax'^y + Axy"^ + 9x|/ + 3y + 6y^ + 
Ay^ + y\ 



TwY^.y) = x"^ + Ax^ + Qx^ + Ax + Ax^y+Axy'^ ^^xy + Ay + Qy"^ + Ay^ + y^ . (46) 





Figure 10: The 4-wheel 



Figure 11: The 4- whirl 



Matroids Qq, Pq, Rq and f/36 

We use (4) to compute the Tutte polynomial of Qq, Pq and Us^, see Figure 12 
for a geometric representation of these matroids. The matroid f/3^6 is uniform, 
from (3) we get 



^c/3.1 



(x, y) = x^ + 3x^ + 6x + 6y + 3y^ + y^ . 
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w^ 



Qe 



Pfi 



U:- 



3,6 



Figure 12: M^K^), W^, Qe, Pg and U; 



3,6 



As f/3,6 is the only relaxation of Pq, from (4) its Tutte polynomial is 

Tpg (x, y) = x^ + 3x^ + 5x + xy + 5t/ + 3y^ + y^. 

Similarly, Pq is the only relaxation of Qq. Thus, from the above equation 
and (4) we obtain 



TQg (x, y) = x^ + 3x^ + 4x + 2xy + ^y + Sy"^ + y^. 



(47) 



It is worth mentioning that Pq is a relaxation of -Rg, see Figure 13. Thus, 
Qq and Rq have the same Tutte polynomial. This is not at all uncommon, 
see [37, 7]. Also, Rq is sparse paving so by (5) its Tutte polynomial is 

TR^{x,y) = (x-l)^ + 6(x- l)2 + 13x- 10 + 2xi/ + 13i/ + 6(i/- 1)^ 
= x^ + 3x^ + 4x + 2xy + 4y + 3y^ + y^. 



Matroids F7, F^ and their duals 

One of the most mentioned matroid in the literature is the Fano matroid, 
Fj, also known as the projective plane PG{2, 2) or the unique Steiner system 
5'(2,3, 7). Being a Steiner triple system implies that the Fano matroid is 
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Figure 13: The matroid Rq 



sparse paving, see [50]. As it has 7 circuit-hyperplanes we obtain from (5) 
that the Tutte polynomial of F^ is 

TF,{x,y) = {x-lf + 7{x-lf + Ux-2l + 7xy + 28y + 2l{y-lf + 
7{y-lf + {y-ir 
= x^ + Ax"^ + 3x + 7xy + 3y + Qy'^ + 3y^ + y'^. (48) 

The matroid F7 is representable over any field of characteristic 2, so the above 
calculation can be checked using the program in [2] with the matrix 

110 1 
10 11 
111 

By using (2) the Tutte polynomial of Fj is 

Tp. {x, y) =x^ + 3x^ + 6x^ + 3x + 7xy + 3y + Ay'^ + y^. (49) 





Figure 14: The Fano matroid and its dual 



The matroids Fj and {Fj )* are the corresponding relaxations of Fj y 
Fy, thus we get that 



3 1 „,4 



Tp- (x, y) = x-^ + Ax^ + 4x + Qxy + Ay + Qy^ + 3y^ + y 



35 





and 



Figure 15: Geometric representation of P-j and Pg 



T(^p~Y (x, y) = x*^ + 3x^ + 6x^ + 4x + 6x|/ + % + ^^ + y^ . 



Matroids P-j,P^ and Qs 

The matroid P7 is a rank-3 sparse paving matroid, thus by using (5) we get 
its Tutte polynomial, 

Tp^{x, y) = X + Ax + 5a: + bxy + by + Qy + 3t/ + y . 

The matrix that represents Py over a field different from GF{2) is 

10 11 
110 1 
a 1 a — 1 

with a ^ {0, 1}. Taking a = 2 we have a representation of P7 over ^^(3), 
see [37]. Thus, the above calculation of the Tutte polynomial can be checked 
using the computer program in [2] . 

The matroid Pg is also sparse paving and its representation over G'P(3) 
is as below, 

011-1" 
10 1 1 
110 1 
-111 

Thus, its Tutte polynomial can be computed using either (5) or the program 
in [2], and you get as a result the following polynomial. 



Tps{x,y) 



X 



4:X^ + lOx^ + lOx + lOxy + lOy + lOy^ + 4y^ + y\ 
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The rank-3 ternary Dowling geometry Qs has geometric representation 
shown in Figure 16. 




Figure 16: Geometric representation of Q? 



The matroid Q^ is representable over a field F if and only if the char- 
acteristic of F is different from 2. A matrix that represents Qs over GF(2>) 
is 

12 3| 4567 89 

"1 0| 1 1 1 1 0' 

1 0| 1 -1 11 

1| 1-1-11 

Thus, we can compute its Tutte polynomial using [2] and obtain 

Tq3 (x, y)=x^ + 6x^ + 8x + 3xy^ + lOxy + 8y + 12y'^ + lOy^ + 6y^ + 3y'' + y^. 

Note, however, that this matroid is paving, so we could have used Proposi- 
tion 2.8 to get the same result. 



Matroid Wl 

The matroid W^l is obtained from W^ by adding an element in parallel and so, 
the matroid is not paving. In this case, we can use Definition 2.3 to compute 
the Tutte polynomial. In Figure 17 we show how we are using deletion and 
contraction to find matroids where the Tutte polynomial is either known or 
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i (7l,3 ffi (/o,2 



Figure 17: Deletion and contraction reduction for W\ 

easy to compute. 

^H/3 (a^, y) = Tw3{x, y) +y- (Tu^A^, y)) + Tu^^^.f^UoA^^ v) 

= {x^ + 3x'^ + 3xy + 3x + 3y + 3y'^ + y^) + y ■ {x"^ + 2x + 2y + y^) 

+ {y'' + y^ + xy^) 
= x^ + 3x'^ + 3a; + x'^y + 5xy + xy"^ + 3y + 5y'^ + 3y^ + y^. 

Matroids ^^(3,2), ^^(3,2)', /?§, Qs, Fs and Lg 

The second affine plane that we consider here is AG{3, 2), Figures 18 and 19 
show two ways of representing the matroid. For example, in Figure 19 the 
planes with 4 points are the six faces of the cube, the six diagonal planes 
{1,2, 7, 8}, {2, 3, 5, 8}, {3, 4, 5, 6}, {1,4, 6, 7}, {1, 3, 5, 7} and {2,4,6,8}, plus 
the two twisted planes {1, 3, 6, 8} and {2, 4, 5, 7}. Each of these planes is a 
circuit-hyperplane. It is not difficult to check that ^^(3, 2) is isomorphic 
to the unique Steiner system S'(3,4, 8), thus, as explained for ylG'(2,3), it is 
sparse paving and its Tutte polynomial is 



T, 



AG(3,2) 



[x, y) =x'^ + ix^ + lOx^ + 6a; + Uxy + 6y + lOy^ + 4y^ + y^. (50) 



By using (4), we can compute the Tutte polynomial of AG{3,2)', the 
unique relaxation of ^^(3, 2) that here we obtained by relaxing the twisted 
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Figure 18: The affine plane 
^^(3,2) 
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Figure 19: The cube with 
points 



plane {2, 4, 5, 7}. Thus, the resuhing polynomial is 

TAG{z,2y{x, y)=x^ \ 4x^ + lOx^ + 7x + I'ixy + 7i/ + lOy^ + 4i/^ + y^. (51) 

Now, AG(3, 2)' has two relaxations, R^ and Fg. The matroid i?8 is ob- 
tained by relaxing from AG(3, 2)' the other twisted plane {1, 3, 6, 8}. On the 
other hand, Fg, is obtained from AG'(3, 2)' by relaxing a diagonal plane. The 
geometric representation of Fg is shown in Figure 20. 




Figure 20: The matroid Fg 

The matroid -Rg is representable over any field except for GF{2) while Fg 
is not representable; however, they have the same Tutte polynomial as both 
are relaxations of the same matroid. 

Tr^ (x, y) = Tf^ {x, y) = x^ + 4x^ + lOa;^ + 8x + 12xy + 8y + lOy^ + 4y^ + y^. 
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The unique relaxation of -Rg is Qs, and is obtained by relaxing one of the 
six diagonal planes of Rg. From the previous equation and (4) we get 

Tq^{x, y) = x^ + 4:X^ + lOx^ + 7x + llxy + 7y + lOy^ + Ay^ + y^. 

To end this subsection we consider the matroid Lg which is a rank-4 sparse 
paving matroid with 8 elements and its circuit-hyperplanes are the six faces of 
the cube plus the two twisted planes {1, 8, 3, 6} and {2, 7, 4, 5}, see Figure 19. 
Thus, its Tutte polynomial is 

TL,{x,y) = (a;-l)^ + 8(a;- l)^ + 28(a;- l)^ + 48a;-42 + 8a;|/ + 48|/ + 
28{y-l)' + 8{y-lf + {y-l)'' 
= x^ + 4x^ + 10a;2 + 12x + 8xy + 12y + lOy^ + 4y^ + y'^. 

Matroids S'g,T8 and J 

The matroid Ss has a geometric representation shown in Figure 21 together 
with its representation over GF{2). Its Tutte polynomial can be computed 
using the program in [2] . 

Tss{x,y) =x^+4:X^ + 7x'^+4:X+10xy + 3xy^ + 3x'^y+4y+7y'^+4y^+y^. (52) 

Note that the matroid is self-dual but is not paving as it has a 3-circuit. 
However, we can check the above computation using deletion and contraction. 
When we contract the forth column in the representation we obtain the 
representation of Fj, and we have already computed the Tutte polynomial 
of this matroid. Now, when we delete the same element, we obtain a rank 4 
graphic matroid. The graph is i^2,4 with an edge contracted. If we call this 
graph H, by using (7) we obtain Th = 7/^2 4 ~ 1/^/^2 a- These polynomial can 
be computed either using the general method for complete bipartite graphs 
given in (30) or by using the formula for computing the Tutte polynomial of 
the 2-stretching of the graphs with two and three parallel edges respectively, 
given in (41). In both cases, we get 

Th = x^ + 3x^ + 3x'^ + X + 3xy + 3x'^y + 3xy'^ + y + y'^ + y^. (53) 

By adding (53) and (48) we get (52) 

The matroid Tg is representable over a field F if and only if the characteris- 
tic is 3. We show its geometric representation in Figure 22. A representation 
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Figure 21: On the left hand side we give the geometric representation of the 
matroid 5*8 and on the right hand side the matrix that represents it over 
GF{2). 



of Tg over GF{3) is [/4I J4 — I^], where J4 is the matrix of I's. This matroid 
is self-dual and sparse paving so its Tutte polynomial is 

Tts{x, y) =x'^ + 4:X^ + 10a;2 + 9x + llxy + 9y + lOy'^ + 4y^ + y^. 




Figure 22: The matroid T8 



In Figure 23 we show the geometric representation of the matroid J. 
It is a self-dual matroid that is not paving, as it has a 3-circuit, and its 
representation over GF{3) is 



10 1 
1110 
10 1 
11 
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where the labelhng of the columns correspond to the labelhng of the elements 
in the geometric representation of Figure 23. The Tutte polynomial can be 
computed using [2]. 

Tj{x, y) =x^ + Ax^ + 7x^ + 6x + 3x^y + 3xy^ + 8xy + Qy + 7?/^ + 4y^ + y*. 




Figure 23: The matroid J 

This computation can be checked by, for example, contracting the element 
labelled 2 to obtain the graphic matroid of the graph that is the 2-stretching 
of C4 minus an edge. Now, the matroid J \ 2 is sparse paving and contains 5 
circuit-hyperplanes. By adding the Tutte polynomial of these two matroids 
you get the same result as above. 



The Vamos matroid 

The 8-element rank-4 matroid whose geometric representation is shown in 
Figure 24 is known as the Vdmos matroid and it is usually denoted by Vs. It 
occurs quite frequently in Oxley's book and has many interesting properties, 
for example is not representable over any field and it is sparse paving, see [37]. 
A related matroid Vg"^ has the same ground set and is defined in the same 
way as Vs but with {5, 6, 7, 8} added as a hyperplane. In fact, Vs is obtained 
from V^g^ by relaxing {5,6,7,8}. Notice that both are self-dual. Thus the 
Tutte polynomial of these two matroids are 



Tvs{x,y) 



(x-ir + 



[X 



11 





V. 



(y-if 



+ 5{xy — X ~ y) + 

+8(t/ - If + (y - ir 

x"^ + Ax^ + lOa;^ + 15x + 5xy + 15y + lOy^ + 4y^ + y'^ (54) 
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and 



Ty+ (x, y) = x^ + 4a;^ + lOx^ + Ux + 6xy + Uy + lOy^ + 4y^ + t/^ (55) 




Figure 24: Geometric representation of Vs 



Matroids Rg„ Rg, Riq and R 
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The real affine cube, Rg, is represented over all fields of characteristic other 
than two by the matrix 
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1 


1 


1 


1 


-1 


1 


1 


1 


1 


-1 


1 


1 


1 


1 


-1 



The matroid is sparse paving and its Tutte polynomial is 

Tr^ {x, y) =x* + 4x^ + lOx^ + 8x + 12xy + 8y + lOy^ + Ay^ + y^. (56) 

On the contrary, the ternary Reid geometry, Rg, is not paving but it is 
representable if and only if the characteristic of the field is three. The matrix 
that represents Rg over GF(3) is 



1 


-1 


-1 


-1 1 











1 


-1 1 


-1 



and its Tutte polynomial can be computed using the program in [2]. 



TR,{x,y) 



X 



6x^ + 8x+nxy + 2xy^ + 8y + 13y^ + lOy^ + 6y^ + 3y^ + y^ 
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This was computed before in section 2.3 with the same result. 

The unique 10-element regular matroid that is neither graphic or co- 
graphic, i?io, has the property that any every single-element deletion is iso- 
morphic to M(i^3 3), and every single-element contraction is isomorphic to 
M*(i^3^3), then by using the polynomial in (31) and (2) we obtain 



TR,^Xx,y) 



x^ + 5x^ + 15x^ + 20a;^ + lOx 
lOy + 20y^ + I5y^ + 5y^ + y\ 



15x'^y + 30xy + 15xy'^ + 



(57) 



Another important regular matroid that is neither graphic or cographic 
is Ri2, which has a matrix representation over GF{2) given by 
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1 











1 


1 





1 








1 











1 








1 











1 








1 





1 


1 











1 


1 


1 



As this matroid is not paving we use the program in [2] to compute its Tutte 
polynomial. 

Tr^^ [x, y) = x^ + 6x^ + 19a:^ + 35a;^ + 35x^ + Ux + 2xS + IQx^y + 
53x^y + Ux'^y^ + 5Qxy + 53xy^ + 19xy^ + 2xy^ + 1% + 
35i/2 + 3hy^ + 19i/^ + Qy^ + y^ 



Pappus and non-Pappus matroids 

The geometric representation of the Pappus matroid is shown in Figure 25. 
From the picture it is clear that any two points are in a unique line and 
that each line is a circuit-hyperplane in a rank-3 matroid. We conclude that 
the matroid is sparse paving with A = 9 circuit-hyperplanes and its Tutte 
polynomial is 



TM{x,y) 



(X 



1)^ + 9{x - 1)2 + 27a; - 78 + 9xy + 117y + 126{y - if + 



If + (y 



84{y~lf + 3Q{y-l)^ + 9{y 

x^ + 6x^ + 12a; + 9xy + 12y + 15y^ + lOy^ + 6y'^ + 3y^ + y^ 
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The non-Pappus matroid is a relaxation of the Pappus matroid and its 
geometric representation is shown in Figure 26. From the previous equation 
and (4) we obtained 

Tm{x, y) = x^ + 6x^ + 12a; + 9xy + 12y + 15y^ + lOy^ + 6y^ + 3y^ + y^ 
—xy + X + y 
= x^ + 6x^ + 13a; + 8xy + 13y + 15y^ + lOy^ + 6y^ + 3y^ + y^. 





Figure 25: Geometric repre- 
sentation of the Pappus ma- 
troid 



Figure 26: Geometric repre- 
sentation of the non-Pappus 
matroid 



Matroid non-Desargues 

The non-Desargues has the geometric representation shown in Figure 27. 
It is a rank-3 matroid with 10 elements that is sparse paving with the 9 
circuit-hyperplanes that are the 9 lines in the picture. Its Tutte polynomial 
is 



TM{x,y) 



x^+7x^ + l9x+9xy+l9y+2ly^~ 



-15y^ + my^+6y^+3y^+y\ (58) 



The Steiner systems 5(2,3,13) and 5(5,6,12) 

Finally, we consider two Steiner systems from the Appendix of Oxley's book. 
In general any Steiner system S{v — l,f,n) is sparse paving, see [50]. Thus, 
the Steiner triple system S{2, 3, 13) is a rank-3 sparse paving matroid with 
13 elements and A = 26 circuit-hyperplanes, then its Tutte polynomial is 



T. 



5(2,3,13) 



[x,y) 



x^ + lOa;^ 



29x 
+28y^ + 21y^ + I5y^ 



26xy + 29y + A5y^ + 3Qy^ 



Wy^ + Qy"^ + 3y^ + y'^ . (59) 
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Figure 27: Geometric representation of the non-Desargues matroid 



The Steiner system S'(5, 6, 12) has rank 6, 12 elements and 132 circuit- 
hyperplanes. The Tutte polynomial is, by using equation (5), 

Ts{5,6A2) {x, y) = x^ + 6x^ + 21a;^ + 56x^ + 126x2 + 120a; + 132xy 

+I20y + 126y^ + 56y^ + 2ly^ + 6y^ + y^. (60) 
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